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^ \ Abstract 

The competition between the isovector [T = 1) and isoscalar (T = 0) proton- 
neutron {p — n) correlations in = Z nuclei is investigated by calculating their 
^ , correlation energies with a realistic effective interaction which reproduces observed 

' nuclear properties very well, in a strict shell model treatment. It is shown in the 

■ realistic shell model that the double-differences of binding energies {B{A + pn : 
^ ' T) + B{A) - B{A + p) - B{A + n)) {B{A) being the binding energy) are good 

indicators of the T = 1 and T = p — n correlations. Each of them, however, 

■ originates in plural kinds of correlations with T = 1 or T = 0. 

PACS: 21.10.-k,21.10.Dr,21.60.Cs 



1 Introduction 



X. 

^ , The competition between isovector (T = 1) and isoscalar (T = 0) pairing cor- 

relations has been a matter of renewed concern in nuclear structure studies of 
N ^ Z nuclei [1,2,3]. The competition appears in the near degeneracy of the 
lowest T = 1 and T = states in odd-odd N = Z nuclei. The T = proton- 
neutron {p — n) pairing correlations in N Z nuclei have been studied by the 
two approaches [2,3] with different treatments of the symmetry energy. The 
two conclusions about the importance of the T = p — n pairing correlations 
are in opposition to one another. The T = 1 and T = pairing correlations, 
which are considered to be induced by T = 1 and T = nuclear interactions, 
should be treated consistently on the same footing [3] . The structure of N = Z 
nuclei has been well described by the shell model which does treat T = 1 and 
T = pairing consistently. Large-scale shell model calculations were applied 
to the investigation of the isovector and isoscalar pairing correlations in Refs. 
[4,5,6], where the contributions of T = 1, J = and T = 0, J = 1 interactions 
are compared [4] and the contributions of the quadrupole-quadrupole (QQ) 
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force are also considered [5]. The authors have recently shown that the com- 
petition between the T = 1 and T = pairing correlations are approximately 
explained with the T = 1, J = pairing force (Pq) and a T = p — n force 
i^mp^ below) [7]. In order to understand the competition in detail, it is im- 
portant to evaluate the two types of correlations induced by realistic effective 
interactions. 

In this paper, we investigate the competition between T = 1 and T = p — n 
correlations in the lowest states oi N — Z nuclei using a reahstic effective 
interaction in a strict shell model treatment. The spherical shell model, which 
gives an excellent description of various properties of N ^ Z nuclei (not only 
the binding energies but also other properties), has the advantage that the 
correlation energies of respective interactions are properly calculated. A shell 
model Hamiltonian is composed of the T — and T — 1 interactions. 



H^Hsp + V^=^ + V^=\ (1) 
V^^J2JlJlJlGjT{ab:cd)Alj^TK{ab)AjMTK{cd) (T = 0, 1), (2) 

a<b c<d JM TK 

where Hgp stands for the single-particle energies, ^jmtx('^^) the creation 
operator of a nucleon pair with the spin JM and the isospin TK in the single- 
particle orbits (a, 6), and Gjriab : cd) denotes the interaction matrix elements. 
The so-called realistic effective interactions contain multipole ( J > 0) pairing 
forces of T = and T = 1 in the expression (2). In this sense, the shell model 
with a realistic effective interaction is to deal with all the multipole pairing 
correlations. We investigate the competition between the T = and T — 1 
correlations induced by the T — and T — 1 interactions. 

The realistic effective interactions have the property that the centroid of T = 
diagonal interaction matrix elements GT=o{ab) = J2j{'^J + l)Gjo(a6 : ab) 
I Y.ji'^J + 1) has a roughly constant value, being independent on the orbits 
(ab) [8]. By setting —ko — Sab G^T=o(ofe)/ Z^ab, we obtain the average T — 
p — n force 

= -^°EE4Moo(«&)^JMoo(a&). (3) 

a<hJM 



Let US write residual T = interactions as V^^s"° = — V^~^ and rewrite 
the Hamiltonian as 

H-H,p + v^;' + vZ7' + v^-\ (4) 



The separation of ° in Eq. (4) follows the procedure of Dufour and Zuker 
in Ref. [9], where the Hamiltonian is divided into the monopole and multi- 
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pole parts as H = Hm + Hm- The monopole field V^p ° is exactly written 
as —{k^/2){n^/2{hj;/2 + 1) — T(T + 1)}, where n„ stands for the number of 
valence nucleons and T stands for the total isospin. This equation shows that 
the force V^p"° which yields only energy gain depending on riy and T is com- 
pletely irrelevant to wave-functions, even to whether the nucleus is spherical 
or deformed, and configuration mixing is caused by V^^^^ and V'^^^. The force 
strength extracted from the realistic effective interaction USD [10] is, for 
instance, 2.8 MeV for ^^Na. We know that Vj^p^ makes a large contribution to 
the binding energy and gives a bonus 2.8 MeV to the T = states against the 
T — 1 states in ^^Na. In this paper, we do not separate the T — 1 monopole 
field which is much smaller than the T = one (we shall mention another 
reason for it later). We analyze roles of the three interactions V^^^, Klr^ 
V'^^^ in the competition between the T = and T = 1 correlations. 



2 Indicators of T = and T = 1 p — n correlations 

It is well known that the even-even N = Z nuclei are more deeply bound than 
the T = 1, O"*" states of neighboring nuclei. The difference in mass between 
the nuclei with mass number A = ma and A = ma + 2 (where the unit a 
consists of two protons (2p) and two neutrons (2n) and m is an integer) is a 
good indicator of the a- like 2p — 2n correlations according to Ref. [11]. If the 
difference is large, it shows the a-likc superfluidity of the A = ma system. In 
fact, we see such a signature in the whole region of N ^ Z nuclei. We can 
regard the T = 1 (or T = 0) lowest state of the odd-odd N = Z nuclei as a 
correlated state of the last p — n pair coupled with the a- like superfiuid state of 
A — ma. Based on this picture, we attempted to analyze the single-difference 
of binding energies B{ma+pn) — B{ma), but we were unable to obtain simple 
information about the T = 1 or T = correlations of the last p — n pair. 

The p — n correlation energy in the odd-odd N — Z nuclei with A — ma + pn 
can be evaluated by the double-difference of binding energies [12,13,14] 

- 5%n{A + pn : T) = {B{A+pn : T) + B{A)) 

-{B{A+p) + B{A + n)) (T = 0,1), (5) 

where B{A) stands for the binding energy. It is popular to calculate the odd- 
even mass difference for evaluation of the T = 1 pairing correlations. Using 
the odd-even mass difference AMoe(A-l-n), the correlation energy of the T —1 
neutron pair in the A + 2n (even) nuclei is evaluated by 

-5V„„(A + 2n) = -2AM,e(A + n) 
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= B{A + 2n) + B{A) - 2B{A + n). 



(6) 



In the BCS approximation for the T = 1, J = pairing correlations, the odd- 
A system is described as one neutron quasi-particle on the pairing superfluid 
vacuum as the average state of adjacent even-A systems, and the quasi-particle 
energy Eqp is approximately equal to the gap A„, i.e., AMoe ~ Eqp A„. 
Since the isospin is a good quantum number in N ^ Z nuclei, the ma + 2p, 
ma + 2n and ma + pn nuclei with T = 1 (or ma + p and ma + n nuclei) have 
approximately the same energy if the Coulomb energy is subtracted from the 
binding energy. We have the approximate relation —S'^Vpn{ma+pn : T — 1) 
— 5^T4„(mQ; + 2n). We can suppose that the quantity —5^Vpn{ma+pn : T = 1) 
represents the correlation energy of the last T = 1 p — n pair as —5'^Vnn{fT^(^ + 
2n) does that of the last neutron pair. Similarly, —6^Vpn{ma + pn : T = 0) 
represents the correlation energy of the last p — n pair with T — outside the 
A — ma system. 

Let us illustrate empirical values of —S'^Vpn{T = 0) and — 5^V^„(T = 1) for the 
lowest T = and T = 1 states of odd-odd N = Z nuclei, in Fig. 1. The dif- 
ference 5'^Vpn{T = 0) — S'^Vpn{T = 1) is equal to the energy difference between 
the lowest T — and T = 1 states, from the definition (5). It is considered 
in phenomenology that the quantities — 5^V^,„(T = 0) and —S'^VpniT = 1) are 
indicators of the T = and T = 1 pairing correlations. The odd-even mass 
difference expressed as 5.18^4"^/^ in the mass formula [15] corresponds well 
with 5'^Vpn{T — l)/2, as shown in Fig. 1. If we write the symmetry energy 
plus Wigner energy as aTT{T -|- 1) in a mass formula, we have the relation 
5'^Vpn{T = 0) ~ 3aT/2 (note that if the Wigner energy proportional to T is ne- 
glected, 5'^VpniT = 0) ^ aT/2). The empirical value ar ~ 134.4(1-1.52^-^/3) 
in the mass formula [15] explains the A-dependent value of 5'^Vpn{T = 0) in 
Fig. 1. The empirical values of 5'^Vp,XT = 1) and 6'^Vpn{T = 0) show large 








18 26 34 42 50 58 66 A 



Fig. 1. Empirical values of —6'^Vpn for the lowest T = and T = 1 states of odd-odd 
N = Z nuclei. 
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fluctuations, which indicates shell effects. Basically, however, S'^Vpn{T = 0) is 
larger than 5^l^„(T = 1) in the sd shell nuclei, while the latter is larger than 
the former in the pf shell nuclei. The empirical trend depending on the mass 
A is well explained by using the parameters of the mass formula [15]. The 
trend of S^Vpn{T — 0) is also reproduced by Sar/^ with somewhat different 
parameters of other modern mass formulas which have the symmetry energy 
and Wigner energy in the form aTT{T + 1). This suggests an intimate relation 
between the T = correlations and the symmetry energy. The model calcu- 
lations using the T — 1 pairing force and a T — p — n force on a deformed 
base [3] and on a spherical base [7] gave a microscopic explanation of the com- 
petition between the T = 1 and T = correlations, reproducing the energy 
difference S%n{T = 0) - S^Vpn{T = 1). 

Our subject is to analyze in more detail the p—n correlation energies —5^Vpn{T — 
0) and —6'^Vpn{T — 1) in the shell model calculations with realistic effective 
interactions. Let us illustrate the results obtained with the USD interaction 
[10] for the srf-shell nuclei, in Fig. 2, where contributions from V^^^, V^^^ 
and to —S^Vpn are denoted by the dash, dot and solid lines respectively, 
and —5'^ {Hap) denotes the contribution from Hgp. The double-difference of 
single-particle energies —5^{Hsp) is evaluated by substituting the expectation 
value (Hsp) for B in Eq. (5). The quantity —S'^{Hsp) gives the starting point 
of —S'^Vpn with no interactions in the shell model calculation. In Fig. 2, the 
calculated values of —d^Vpn{T = 0) and —S'^Vpn{T = 1) finely reproduce their 
empirical values. We see mixed contributions from the T = and T — 1 
interactions to the final states with T = or T = 1 and some deviations 
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Fig. 2. The double-differences of binding energies —S^Vpn{ma + pn : T) and 

— (5^V2p2n(f^o+2p2n). The starting point of —S'^Vpn or — 5^V2p2n with no interactions 
is shown by the top fine (—6'^ (Hsp)) in each column. The respective contributions 
from V^p^, V^^7° and V"^^^ are denoted by the dash, dot and solid lines. The open 
diamonds denote the empirical values. 
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coming from the single-particle energy part Hgp. Still, Fig. 2 shows that the 
double-differences of binding energies — 5^l^„(T = 0) and —6'^Vpn{T = 1) are 
good indicators of the T = and T = 1 p — n correlations. It is notable that 
not only V^^^ but also V^^^ contributes to —5'^Vpn{T = 0). We calculated 
separately the contributions from the T — 1, J — interactions and the other 
(T = 1, J > 0) interactions to the T = 1 p — n correlation energy. The cal- 
culated indicator — 5^V^„(T = 1) includes the contributions from the T = 1, 
J > interactions as well as the T = 1, J = interactions. The results show 
that the T = 1, J > interactions contribute greatly to the binding energy 
but reduce the value of —6'^Vpn{T — 1) considerably. If we consider only the 
T = 1, J = pairing force Pq and adjust the force strength so as to explain 
the quantities 6'^Vpn{T = 1) and 6'^Vnn for the last nucleon pair, the binding 
energy will not be reproduced, revealing a contradiction. 

We can evaluate the p — n correlation energy in even-even N — Z nuclei using 
the following double-difference of binding energies [8,14,16,17,18]: 



- SW2p2„{A + 2p2n) = {B{A + 2p2n) + B{A)) 

-{B{A + 2p) + B{A + 2n)). (7) 

When we consider the A = ma systems, the quantity S'^V2p2n is nothing but the 
difference in mass caused by the a-hke correlations [11,19]. The large values 
of 5^V2p2n observed in the even-even N — Z nuclei show the cc-like superfluid- 
ity. The definition of S'^V2p2n based on the picture of the a-like superfluidity 
corresponds well with that of S'^Vnn based on the pairing superfluidity. The 
quantity —S^V2p2n represents the p — n correlation energy between 2p and 2n. 
A quarter of —S'^V2p2n (— ^^^p2n/4) is regarded as the p — n correlation energy 
per one p — n bond. 

In Fig. 2, we show the calculated and empirical values of — 5^V2p2n/4 for the sd- 
shell nuclei ^"^Mg and ^®Si. This flgure indicates that the quantity — 5^V2p2n/4 
evaluates the T = p — n correlation energy as the quantity —6^Vpn{T = 0). 
The value of S'^V2p2n/4: is smaller than that of 5'^Vpn{T = 0), which means that 
the p — n correlations of the last p — n pair with T = in the A — ma + pn 
nuclei are a little stronger than the p — n correlations between 2p and 2n. 
By the way, 5^V2p2n corresponds to the symmetry energy of the A + 2p or 
A + 2n system with T = 1 in the framework of the mass formula [20] . Thus, 
the indicators S'^Vpn{T = 0) and S'^V2p2n represent the T — p — n correlation 
energy in the shell model, while they are attributed to the symmetry energy 
in the mass formulas. If one excludes the symmetry energy from the binding 
energy as done in Ref. [2], one inevitably misses the contribution from the 
T — p — n interactions to the total energy and hence is likely to miss the 
signature of the T = correlations. 
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3 Further discussions of T = and T — 1 correlations 



It is instructive for our study of the T — 1 and T — correlations to see the ex- 
pectation values of V^^'^, Klr° V"^^^ in addition to the double-differences 
of binding energies. We calculated their expectation values in the lowest states 
with T = and T = 1. We call the expectation values "interaction energies" 
and write them such as {V'^^^). The interaction energies obtained by using 
the USD interaction for the N — Z sd-shell nuclei from ^°Ne to ^^S are shown 
with three different lines in Fig. 3. This figure shows the indivisible coopera- 
tion of the T = 1 correlations and the T = p — n correlations, and indicates 
a significant role of the average T = p — n force V^^^ in the realistic ef- 
fective interaction. The larger the valence-nucleon number is, the greater is 
the V^p^ contribution to the binding energy. Figure 3 shows that the T — 
p — n interactions including V^p^ make a larger contribution to the binding 
energy than the T = 1 interactions. However, if wc neglect V^f^p^ which does 
not affect the structure, the T = 1 interactions are more important than the 
T — residual interactions. 



Dufour and Zuker [9] indicated that the realistic effective interactions are 
represented approximately by the monopole field V„ip, the Pq force and the 
QQ force. In fact, it was shown in Ref. [21] that the interaction Vmp + Po + QQ 
with the T — 1, J — 2 pairing force (P2) reproduces various properties of 
the /7/2-subshell nuclei as comparably as the realistic effective interactions. 
The extended P + QQ interaction describes well also the s(i-shell nuclei in 
our calculations (unpublished). These results showed that the essential part 
of the monopole field is the T = one (V^p^) and the T = 1 monopole field 
can be neglected, although some monopole corrections (AV^""° and AV^z"^) 
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Fig. 3. Interaction energies of V^p °, V^J^g ° and ^ in the sd-shell nuclei. Their 
sum is the total interaction energy. 
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depending on respective orbits remain. This is another reason why we do not 
extract the T = 1 monopole field in Eq. (4). The reahstic effective interactions 
are well approximated by 

v^;' + Av^;' + Av^;' + P0 + P2 + QQ. (8) 

We have a useful insight into the T = and T = 1 correlations in this model. 
The T = correlations are induced by T = re-coupling terms of the QQ 
force and AV^p^ except for V^p^, while the T = 1 correlations are induced by 
the monopole and quadrupole pairing forces (Pq and P2), and also by T = 1 
re-coupling terms of the QQ force and AV^p^. 

As is well-known, the T = states are the ground states in the odd-odd N = Z 
nuclei except for ^^Cl in the sd shell. It is curious that the interaction energy 
gain of the T — 1, O''' state is larger than that of the lowest T = state for 
2^A1 and in Fig. 3. In the shell model, the T = states of ^^A\ and 
suffer less energy loss due to the single-particle energy part (Hsp) and hence 
become the ground states. This happens because the T = 1 correlations are 
stronger than the T = residual p — n ones, thus making more nucleons jump 
up to the ds/2 orbit. If we lower the ds/2 orbit, the T — 1 state becomes lowest 
in ^^Al. The shell structure thus affects the competition between the T — 
and T — 1 states for the ground-state position in the odd-odd N — Z nuclei. 



4 Conclusion 

In conclusion, we have investigated the competition of the T = and T = 1 p— 
n correlations in an exact treatment of the shell model with a realistic effective 
interaction. By dividing the T — p—n interactions into the average monopole 
field and the residual interactions, we evaluated interaction energies of them 
and the T = 1 interactions. The calculations show that the cooperation of 
the average T = p — n force, T = residual p — n interactions and T = 1 
interactions produce the large binding energies of the N = Z nuclei. It is shown 
that the double-difference of binding energies —S^Vpn{T — 0) or —5'^Vpn{T = 1) 
defined by Eq. (5) can be used to evaluate the correlation energy of the last 
p — n pair with T = or T = 1 in the odd-odd N = Z nuclei. The realistic 
shell model calculations, however, show that the indicator — (5^V^„(T = 0) or 
—S'^Vpn{T = 1) does not originate in a single kind of pairing correlations, but 
contains the contributions from plural kinds of correlations with T = or 
T= 1. 

This work was stimulated by the workshop on "isoscalar and isovector pair- 
ing" organized by Prof. R. Wyss in Stockholm in May 2003. The authors are 
grateful to Dr. D. M. Brink for his help in improving the manuscript. 
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